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The mean force exerted upon a perfect mirror moving in 
vacuum in a two dimensional spacetime has the same expres- 
sion as the radiation reaction force computed in classical elec- 
tron theory. It follows that unacceptable runaway solutions 
are predicted. We show that this instability problem does 
not appear when partially transmitting mirrors are studied. 
The mechanical impedance describing the mirror coupled to 
vacuum radiation pressure is computed explicitly; recoil is ne- 
glected. It is found to be a passive function, so that stability 
is ensured. This is connected to the fact that no energy can 
be extracted from the vacuum state. 



INTRODUCTION 

A mirror moving in vacuum experiences a mean force 
which is associated through linear response theory with 
the quantum fluctuations of the vacuum stress tensor jlj . 

This motional force, denoted hereafter F m , has been 
computed by Fulling and Davies ^| for a perfectly re- 
flecting mirror scattering a scalar vacuum field in a two 
dimensional spacetime. A linear approximation in the 
mirror's displacement leads from their result to a force 
proportional to the third derivative of the time dependent 
mirror's position q 



F m {t) 
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(1) 



As required by Lorentz invariance of the vacuum state, 
this force vanishes for a uniform velocity. It also vanishes 
for a uniform acceleration so that the inertial mass of the 
mirror is not modified. 

The motional force modifies the mechanical response 
of the mirror to an applied force F a since the equation of 
motion can be written 



kq(t)+mq"(t) = F a {t) + F m {t) 



(2) 



where m is the mass of the mirror, loq is the eigenfre- 
quency of the suspension system (lu = in the limit- 
ing case of a free mass) and F m is the motional force. 



Consequently, it has to be taken into account in a com- 
plete analysis of the sensitivity of interferometric mea- 
surements. This seems particularly important for inter- 
ferometric detection of gravitational waves || . 

The force (1) has the same expression, with a mod- 
ification of the numerical factor, as the radiation reac- 
tion force computed in classical electron theory [Q. It 
is well known that severe problems are associated with 
the resulting equation of motion (2): a historical account 
is given in ref. more recent references can be found 
in p|. First, it turns out that 'runaway solutions' ex- 
ist, which correspond to an exponential increase of the 
velocity in the absence of applied force. These runaway 
solutions may be prohibited by a special choice of bound- 
ary conditions, but a 'preacceleration' phenomenon then 
occurs, with an acceleration of the mass before the force 
is applied 

Clearly, the equation of motion (2) violates causality 
for a perfect mirror when the motional force is given by 
equation (1). We show in this letter that this problem is 
solved in a simple manner by considering partially trans- 
mitting mirrors rather than perfect ones. We will assume 
that the mirror is transparent at frequencies higher than 
a reflection cutoff cue and that this cutoff obeys the fol- 
lowing inequality 



hujc •C mc 2 



(3) 



This condition, obeyed by realistic mirrors (to is a macro- 
scopic mass), implies that the recoil effect for reflected 
fields can be ignored in the derivation of the motional 
force []. 

Using a previously derived expression for the motional 
force associated with a partially transmitting mirror [Q, 
we will demonstrate that the mirror coupled to the vac- 
uum radiation pressure is a stable system (runaway so- 
lutions disappear). More precisely, we will show that the 
mirror obeys 'passivity' Q when some inequality which 
is a consequence of condition (3) is satisfied. Passivity 
is closely connected to energy considerations: the mov- 
ing mirror radiates energy into vacuum and is therefore 
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The perfect mirror corresponds to u>c — * oo and does not 
obey condition (3). It is therefore necessary to take the re- 
coil into account, which was not done in the derivation of the 
force (1). This probably explains the inconsistencies associ- 
ated with expression (1) of the force. 
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damped, while the vacuum cannot provide energy for ex- 
citing the mirror's motion. Then, stability follows and 
problems associated with runaway solutions are elimi- 
nated. 



THE SUSCEPTIBILITY FUNCTION 
ASSOCIATED WITH THE MOTIONAL FORCE 

At the limit of small displacements, the motional force 
is a function of q easily defined in the frequency domain 
in terms of a linear susceptibility \ 



F m (t) = Jdt'x(t-t')q(t') 



(4a) 
(4b) 



Any function / is written in the time domain or in the 
frequency domain according to the rule 



/(*) = / £/Me" 



The motional force (1) computed for a perfect mirror 
corresponds to the susceptibility 



XM 
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(5) 



where r is a time constant, characteristic of the coupling 
of the mirror of mass m to the vacuum radiation pressure. 
The corresponding time constant for the electron is of 
the order of 10~ 22 s. It is much smaller for a macroscopic 
mirror. 

In a scattering approach 0, a partially transmitting 
mirror is described by a reflectivity function r[oj\ and 
a transmitivity function s[w], which obey unitarity and 
causality conditions. The mirror is supposed to become 
transparent at the high frequency limit (r — > when 
w — > oo). The perfect mirror is then replaced by the 
more realistic model of a causal mirror having a perfect 
reflection below a reflection cutoff u>c, but transparent 
at higher frequencies. 

This approach has provided us with an expression for 
the mean Casimir force between two motionless mirrors 
, which is free from the divergences usually associated 
with the infiniteness of the vacuum energy, as well as a 
causal expression for the motional Casimir force between 
moving mirrors 

When condition (3) is obeyed, it is possible to ignore 
the recoil effect and the motional force can be written as 



= imTU 3 r[Lj] 



(6a) 



The motional force for a mirror in vacuum is given by equa- 
tions (33) and (11) of ref. § 



where 

w 3 I>] = / 3du/ [uj - J )uj'a\u - <J , w'] (6b) 
Jo 

a\u,uj'] = 1 - s[w]s^ + r[w]r[w'] (6c) 

The function T[lj] is regular at the low frequency limit 

T[lj] -> T[0] = r[0] 2 for lu -> 

When r[0] = —1, the susceptibility (5) is recovered at 
low frequencies whatever the reflectivity may be at high 
frequencies. At high frequencies, T[u>] behaves as a cutoff 
function. Its behaviour will be discussed in more detail 
later on. 

A simple example is provided by the following model 
fulfilling unitarity, causality and transparency conditions 

r M = 73!^ *[w] = l+r[w] (7a) 

One obtains in this case 

iu> u> 2 
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for u> <C 



(7b) 



POSITIVITY OF THE DISSIPATIVE 
FUNCTIONS 

Coming back to the general case, we introduce the 
real and imaginary parts of the functions written in the 
frequency domain, which are associated with dissipation 
and dispersion 

xH = XrM + ixi[u\ r H = r fiH + iTi[u)\ 
= TOTa; 3 r fl [w] = -m™ 3 r i[uj] 

Xr and Tr are even functions of ui while xi an d Tj are 
odd functions of u (see eqs 4, 6). The dissipative part Xi 
of the susceptibility, which was denoted £f f i n ref. Q , is 
the commutator of the force operator. 

Using the unitarity of the 5— matrix associated with 
the mirror Q , one shows that 

2Rea[w,u'] = \a[uj,uj']\ 2 + \P[uj,lo'}\ 2 
P[u,lo'} = s[u]r[u'] - r[u)]a[u'] 

The following positivity property follows 

r_R[w] > for u) real 

Xi[u] > for u real > (8) 

This property corresponds to the dissipative character 
(see for example ref. Q) of the motional perturbation 
in the vacuum and will play an important role in the 
following. 
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DISPERSION RELATIONS 

Causality |l is associated with the fact that the sus- 
ceptibility function, written in the time domain (see eq 
4), vanishes for negative times 



x(t) = 



for t < 



In the frequency domain, a causal function is analytic and 
regular in the half plane Imcj > 0. For the susceptibility 
given by equations (6), these properties follow from 
the causality of the scattering coefficients r and s. 

The real and imaginary parts of a causal function are 
related through dispersion relations, which take the sim- 
ple form of Kramers-Kronig relations when the function 
decreases sufficiently rapidly at high frequencies || . But 
this is not the case for the susceptibility function Xi so 
that it is necessary to write dispersion relations with sub- 
tractions [Tc| ]. 

From equations (6), one knows that the three quanti- 
ties x[0], x'[0] and x"[0] vanish. It is therefore tempting 
to write the dispersion relation with three subtractions 
at lo = which is just the Kramers-Kronig relation for 
the function T 



ry 



dto' T r [lj'] 



LJ — It 



(9) 



Assuming that T[uj] decreases sufficiently rapidly at high 
frequencies so that it is a square integrable function, the 
relation (9) is actually equivalent to causality for the sus- 
ceptibility function. This follows from the Titchmarsh 
theorem |£o). 

As T R is positive, one deduces that T[ui] behaves as ^ 
at high frequencies 



UJ C 



for 



(10) 



where the parameter loc can be considered as the defini- 
tion of the reflection cutoff (we suppose that the integral 
is finite) 



f dui x/H 

7T / IT mTUJ a 



(11) 



For example, the model of equations (7) corresponds to 
luc = 3f2. 

The behaviour of the susceptibility is described by a 
high frequency mass fi 



/i = mucT 



for uj — > oo 

67TC 2 



(12) 



The induced mass fi is quite similar to the 'electromag- 
netic mass' of electron theory j|. It is positive (see eq 8) 
and much smaller than m for realistic macroscopic mir- 
rors (see eq 3). 



The dispersion relation (9) can be written equivalently 
in the time domain 

r(t) = 26{t)T R (t) 

It is thus clear that the susceptibility function does not 
obey the same simple dispersion relation; one can indeed 
write 

X (t)=2mr(e{t)T%(t)+S(t)TU0) 

+ S'(t)T' R (0)+S"(t)T R (0)) 

It follows that the expression (4a) of the motional force 
cannot be written simply in terms of the force commuta- 
tor xi (t) since quasistatic motions must be subtracted up 
to the second order in a Mac-Laurin expansion. However, 
using the parity property of the function T R , one recovers 
the consistency condition of linear response theory ^ 

x(t)-x(-t) = 2mrT' R \t)=2i Xl (t) 



MECHANICAL IMPEDANCE OF THE MIRROR 

When the coupling of the mirror to vacuum radiation 
pressure is taken into account, the equation of motion (2) 
is easily solved in the frequency domain. 

The mechanical impedance Z gives the frequency com- 
ponents of the force F a as a function of the velocity v 

v(t) = q (t) v[uj] = — iu>q[uj] 
F a [uj] - Z[l>]v[l>] 
—iuiZ[uj] = k — muj 2 — %[cj] (13) 

The two real functions Z R and Zi are respectively the 
dissipative and dispersive parts of the impedance 

Z[uj] = Z r [uj] + iZj[io\ 
Z R [iv] — mTUJ 2 T r [lu] 



Zj[uj] = muj + mTui 2 T 



(14) 



Equivalently, we can write the velocity v as a function 
of the applied force F a and of the mechanical admittance 
Y 



i[w]=Y[w]F a [u] Y[uj} = 



1 



Z[w] 



3 The dispersion relation written in refs. Q and || did not 
account for subtractions and has therefore to be amended. 
However, this does not affect the expressions for the linear 
susceptibility, the force commutator and the force fluctua- 
tions which were independently computed. Also, the conclu- 
sion that these quantities are consistent with linear response 
theory remains valid. 
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STABILITY 

From causality, we know that the motional force is a 
retarded function of the mirror's position, which implies 
that the impedance function Z[uS\ is analytic and has no 
poles in the half plane Imw > 0. 

The analytic properties of the admittance function 
Y[oj] require a closer examination since its expression in 
terms of the susceptibility xM has the form of a 'closed 
loop gain', so that a self oscillation could occur. Such 
a self oscillation would be associated with the presence 
of a pole for Y, that is a zero for Z, in the half plane 
Imw > 0. 

But, if the suspended mirror coupled to the vacuum is 
a stable system, the self oscillation cannot take place. It 
follows that the admittance function Y[u>] shall also be 
analytic and regular in the half plane ImcJ > and that 
the impedance function Z[u>] should have no zero in this 
domain. 

Using the susceptibility function (5) of a perfect mirror 
(with neglected recoil) and assuming that the mirror is 
free (k — 0), one checks that the function Y has a pole 
in the half plane Imut > 0, at ui — -, in this case. This is 
similar to the well known instability problem of electron 
theory g. 

This result is modified when a partially transmitting 
mirror is studied. Indeed, the pole associated with run- 
away solutions corresponds to a frequency well above the 
reflection cutoff (see eq 3). Due to the asymptotic be- 
haviour of the cutoff function T (see eq 10), the pole 
which was located at ui — ^ for a perfect mirror is shifted 
to oj « — iojc so that the instability is suppressed. 



with fj, defined by equations (12) and (11) and the spec- 
tral function $ given by 



d*(p) 



ZR.[p\dp 



One knows that Zr is positive at all frequencies (see eqs 
8 and 14), that the function $ is bounded (this follows 
from the finiteness of p) and that the coefficients k and 
m — fi are positive (when p, < m). Then, equation (16) 
is actually the spectral representation for a passive func- 
tion, that is a function obeying the property |6[ 



ReZ{p} > for Rep > 



(17) 



Causality and stability follow from passivity |l0| , so that 
the problems associated with runaway solutions are elim- 
inated. 

It has to be emphasized that the passivity condition 
(17) is more stringent than the already known positiv- 
ity property (8): as an example, the expressions corre- 
sponding to the perfect mirror obey condition (8) but not 
condition (17). 

It is also worth to stress that the passivity property is 
obtained only when the mass m of the mirror is greater 
than the induced mass p which is similar to an electro- 
magnetic mass. As far as electron theory is concerned, 
the same conclusion is reached by Dekker Jllfl , In other 
words, passivity is obeyed by the impedance Z but not 
by the contribution of the motional susceptibility taken 
separately in contrast with the usual models of Brow- 
nian motion where the motional contribution itself is a 
passive function PJl2| . 



PASSIVITY 



ENERGY CONSIDERATIONS 



We demonstrate now in a more rigorous manner that 
the instability is suppressed for any reflectivity function, 
provided that the following inequality is obeyed 



p < m 



(15) 



In particular, this will be the case as soon as the cutoff 
fulfills condition (3). 

The demonstration relies upon the dispersion relation 
(9) which reads for the impedance function, for Imcj > 



Z[ui] = oj' 



2 f Aoj' Zr[u>'] 



ik 



7r oj' 2 (e — ioj + ioj') OJ 
Using Laplace transforms instead of Fourier transforms 

f{p} = f[ip] Rep > 
this relation can be transformed to 

z{ P }= fm{py—^- + -+p{m-^) (16) 
J p-ip p 



Passivity is equivalent to the proposition that a pos- 
itive amount of energy is transfered from the reservoir 
(source of the applied force F a ) to the moving mirror |j] 



W a (t) = / At' F a {t')v(t') > 



(18) 



One immediately deduces from the equation of motion 
(2) that 

W a (t) = AE(t) + W m (t) 
AE{t) = E(t) - E(-oo) 

E{t) = l -kq{tf + l -mv{tf 



W m (t) 



At' F m {t')v{t') 



4 At low frequencies ~ mr6j 2 r[0], so that this function 
does not obey the passivity property (17) in the vicinity of 
uj = 0. 
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The energy W a provided by the reservoir is the sum of an 
energy AE stored in the mirror's motion and of a part 
W m radiated in the vacuum state. 

When we consider a situation where the force F a is 
applied only during a limited time interval Q and where 
the initial and final states of the mirror are the same 
(AE(oo) = 0), we conclude from condition (18) that the 
net radiated energy is positive 

W a (oo) = W m (oo) > 

Thus, passivity appears to be connected to the fact that 
vacuum can damp the mirror's motion but cannot excite 
it. This is the ultimate reason why stability follows from 
passivity: the vacuum cannot provide energy for sustain- 
ing runaway solutions. 



[8] Jaekel M.T. and Reynaud S., J. Physique I 2 149 (1992). 
[9] Landau L.D. and Lifshitz E.M., Course of Theoretical 
Physics: Statistical Physics (Pergamon, 1958) ch.12. 

[10] Nussenzveig H.M., Causality and dispersion relations 
(Academic Press, 1972). 

[11] Dekker H., Phys. Lett. A107 255 (1985). 

[12] Ford G.W., Kac M. and Mazur P., J. Math. Phys. 6 504 
(1965); Haken H., Rev. Mod. Phys. 47 67 (1975); Dekker 
H., Phys. Rep. 80 1 (1981); Caldeira A.O. and Leggett 
A.J., Ann. Phys. 149 374 (1983); Ford G.W., Lewis J.T. 
and O'Connell R.F., Phys. Rev. Lett. 55 2273 (1985). 



CONCLUSION 



The expression (1) of the damping force for a perfect 
mirror in vacuum is associated with unacceptable run- 
away solutions. We have shown that this problem does 
not take place for a realistic mirror with a reflectivity 
cutoff obeying condition (3). 

The impedance function can be computed explicitly, 
neglecting the recoil effect, and is found to be a passive 
function, so that stability is ensured. This is connected 
to the fact that no energy can be extracted from the 
vacuum state. 

It is tempting to consider that passivity is actually a 
general property of the impedance function for a scatterer 
in vacuum. Then, a consistent treatment of a perfect 
mirror, taking into account the recoil effect for reflected 
fields, would necessarily provide us with a passive func- 
tion. This would also be true for a satisfactory treatment 
of the electron. Attempts in this direction, unsuccessful 
up to now, are discussed by Rueda Q . 
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